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ABSTRACT 

We present photometric and spectroscopic observations of two luminous blue 
variable (LBV) stars in two extremely metal-deficient blue compact dwarf (BCD) 
galaxies, DDO 68 with 12+logO/H = 7.15 and PHL 293B with 12+logO/H = 
7.72. These two BCDs are the lowest-metallicity galaxies where LBV stars have 
been detected, allowing to study the LBV phenomenon in the extremely low 
metallicity regime, and shedding light of the evolution of the first generation 
of massive stars born from primordial gas. We find that the strong outburst 
of the LBV star in DDO 68 occurred sometime between February 2007 and 
January 2008. We have compared the properties of the broad line emission in 
low-met allicity LBVs with those in higher metallicity LBVs. We find that, for 
the LBV star in DDO 68, broad emission with a P Cygni profile is seen in both 
H and He I emission lines. On the other hand, for the LBV star in PHL 293B, 
P Cygni profiles are detected only in H lines. For both LBVs, no heavy element 
emission line such as Fe II was detected. The Ha luminosities of LBV stars in 
both galaxies are comparable to the one obtained for the LBV star in NGC 2363 
(Mrk 71) which has a higher metallicity 12+logO/H = 7.89. On the other hand, 
the terminal velocities of the stellar winds in both low-metallicity LBVs are high, 
~ 800 km s _1 , a factor of ~ 4 higher than the terminal velocities of the winds 
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in high-metallicity LBVs. This suggests that stellar winds at low metallicity are 
driven by a different mechanism than the one operating in high-metallicity winds. 

Subject headings: galaxies: abundances — galaxies: irregular — galaxies: ISM 
— H II regions — stars: winds, outflows 



INTRODUCTION 



The most massive stars evolve on very short timescales of a few million years compared 
to the very long time scales of billion years of solar-type stars. Thus, important evolutionary 
events in the life of very massive stars such a s core-collapse supernovae or giant eruptions of 
luminous blue variables (LBVs; IConti Ill984l ) are hard to catch and study. To maximize the 
probability of catching such an event, it is best to observe objects known to contain many 
massive stars. Blue compact dwarf (BCD) galaxies are ideal laboratories for such studies as 
they are undergoing intense bursts of star formation a nd their star-forming re gions are known 
to harbor up to tens of thousands of massive stars (IThuan fc Izotov 1 120051 ) . Furthermore, 
mass loss is a critical factor in determining the evolution of a massive star. The mass 
loss rate depends in turn on the metallicity of the gas, if driven by line opacity. BCDs 
are also excellent laboratories to test metallicity effects on mass loss as they are metal- 
deficient, their meta llicity ranging from about half to a few percent the Sun's metallicity 
( jlzotov et al. 1 120071 ) . Thus BCDs are the only galaxies in the local Universe where low- 
met allicity massive stars do exist. They are unique objects for studying processes related 
to massive star formation and evolution and their interaction with the interstellar medium, 
and for testing models of massive stellar evolution. They are the best local approximations 
to primordial galaxies. 

BCDs have been used to investigate mass loss through observational signatures of stellar 
winds from massive stars. The prevailing belief is that the efficiencies of stellar winds in 
massive stars are significantly reduced at low metallicities. However, this conclusion is based 
mainly on theoretical extrapolations and observations of relatively high-metallicity massive 
stars in our Galaxy (12 + log O/H ~ 8.7), the Large Magellanic Cloud (12 + log O/H ~ 
8.4) and the Small Magellanic Cloud (12 + log O/H ~ 8.1) and some other nearby high- 
metallicity galaxies. There is growing evidence suggesting that massive stars do have stellar 
winds, even i n the most metal- defi cient BCDs known, with ox ygen abundances 12 + log O/H 
< 7.6. Thus. llzotov et al. I Jl997ah and lLegrand et al. I Jl997h discovered a Wolf-Rayet (WR) 
stellar population in I Zw 18, the most metal- deficient emission-line galax y known at the 
time, with oxygen abundance 12+logO/H = 7.17 ± 0.01 (jlzotov et al. Ill999l ). iGuseva et al 



( I2000l ) have shown that the observed ratio of the number of WR stars to the number of O 
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stars in I Zwl8 cannot be reproduced by theoretical population synthesis models based on 
non-rotating stellar evolution models and observational prop erties of high-metallici t y WR 
stars. To achieve agreement between observations and theory, ICrowther fc Hadfield I (120061 ) 
have suggested that WR stars in I Zw 18 have lower emission-line luminosities than their 
high-metallicity counterparts. Furthermore, iThuan fc Izotov I (119971 ) have observed P Cygni 
profiles for the Si IV A1394, A1403 absorption lines in the spectrum of the emission-line galaxy 
SBS 0335-052E, one of t he most metal-deficient BCD known, with 12 + log O/H = 7.30 ± 
0.01 (jlzotov et al. 1 120071 ). It thus appears that extremely low-metallicity massive stars do 
possess stellar winds, although the properties of these stellar winds may significantly differ 
from those of their high-metallicity counterparts. 

To find more objects with observational signatures of massive star activity, we have 
recently assembled a large sample of about 40 emission-line dwarf galaxies which exhibit 
broa d components in th eir strong emission lines, mainly in H/3, [O ill] AA 4959, 5007, and 
Ha (llzotov et al. 1 120070. Except for fou r objects which appear to contain intermediate- 
mass black holes (jlzotov fc Thuan 1 120081 ) . the broad emission of all the other objects can 
be attributed to some evolutionary stage of massive stars and to their interaction with the 
circumstellar and interstellar medium: WR stars, supernovae, superbubbles or LBV stars. 
In particular, our attention was drawn to the peculiar spectrum of the BCD PHL 293B = 
SDSS J2230-0006, with an oxygen abundance 12 + log O/H = 7.66 ± 0.04, which shows 
broad hydrogen emission lines with P Cygni profiles, spectral featu res that are characteristic 



of a LBV star. Coincidentally, several months after our finding, iPustilnik et al. I (l2008bl ) 
discovered another bright LBV in the BCD DDO 68. This BCD, with an o xygen abundance 
12 + log O/H = 7.14 ± 0.03 is even more metal-deficient than PHL 293B Jlzotov & Thuan 



20071 ). It is the second most metal-deficient star- forming galaxy known, after SBS 0335- 
052W. 

The study of these newly discovered LBV stars in the two very metal- deficient galaxies 
DDO 68 and PHL 293B is the focus of this paper. Giant eruptions of LBV stars, with a 
brightening greater than 3 mag, a re exceedingly rare. In the Galaxy, although there are some 



35 confirmed or candidate LBVs (IClark et al. 1120051 ). only two, rj Carinae and P Cygni, are 



known to have undergone such giant eruptions over the last 5 centuries, rj Carinae between 
1837 and 1860, and P Cygni, twice in the 17th century. P Cygni was discovered in 1600 as a 
naked eye star, and remained bright for many years, before fading and re-appearing in 1655, 
and then fading again. Typically, LBVs exhibit day-to-day microvariations in brightness 
of 0.1-0.2 mag, and "normal" irregular variations of 1-2 mag on a timescale of years, in 



B or early A supergiant at visual maximum 



1994; Leitherer et al. 1994: Nota et al 



(e.g. 


Conti 


1984 




Humphreys & Davidson 


Drissen et al. 


1997. 


2001; 


Crowther 


2004; 
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Petit et al. 1120061 ; IWalborn et al. 1 120081 ) - There ha ve been also LBVs o bserved in the Large 



Magellanic Cloud and other Local Group galaxies (IMassey et al. 1 120071 ) 



The LBV phase is believed to represent a critical transition in the late evolution of 
all stars with initial masses greater than about 50 M & . During this phase, the stars lose 
sufficiently large amounts of mass in recurrent explosive events, to go from the stage of 
O stars burning hydrogen on the main-sequence to that of core-helium burning classical 
WR stars. In the HR diagram , LBVs lie just to the left of the Humphreys-Davidson limit 
( [Humphreys fc Davidson 1119941 ) . beyond which no stars are observed. They appear to define 
the locus of an instability that prevents further redward evolution. 

The understanding of the physical processes operating during the LBV phase is thus 
crucial for modeling massive star evolution. However, the exact mechanism giving rise to a 
LBV outburst is still unknown, even if many LBVs at solar metallicity have been studied. 
The situation is worse at low metallicity because the only relatively low-metallicity LBV star 
that has been observationally investigated in detail is the star VI in the H H region NGC 2363 



(Mrk 71) of the dwarf comet ary galaxy NGC 2366 (IDrissen et al. 



20061), with 12+log O/H = 7.89 flGonzalez-Dekado et al. 



1994 



1997. 



2001 



Petit et al. 



Izotov et al. Ill997bh . High 



quality spectroscopic studies of LBV stars in two more extremely metal-deficient galaxies 
should help us to better understand how the LBV phase depends on metallicity. 

We present here new spectroscopic observations of the LBVs in DDO 68 and PHL 293B 
and derive their emission line parameters. We also discuss new photometric observations 
of DDO 68 to better constrain the epoch of the LBV outburst in it. The observations are 
described in §2. We discuss in §3 the photometric variations of the two LBV stars. We derive 
the element abundances in the two host galaxies and show how the main properties of the 
stellar winds in low-metallicity LBVs differ significantly from those in their high-metallicity 
counterparts. Our conclusions are summarized in §4. 



2. OBSERVATIONS 



2.1. Imaging 



New deep images of DDO 68 were obtained with the SDSS g and % filters ( jFukugita et al. 
1996 ) using the 2.1m telescope at Kitt Peak National Observatoryfj] on the night 9 Febru- 



1 Kitt Peak National Observatory is operated by the Association of Universities for Research in Astronomy 
(AURA) under cooperative agreement with the National Science Foundation. 
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ary 2007. The details of the observations are shown in Table [TJ These observations were 
compared to the SDSS images of DDO 68 in the same filters obtained on 16 April 2004. 
The details of these observations are also shown in Table [TJ To process the 2.1m KPNO 
images we have used the standard data reduction procedures in IRAF§, which include bias 
subtraction and flat-field correction. For photometric calibration, we have used bright stars 
in the field of SDSS images. The resulting g and g — i 2.1m KPNO images of DDO 68 are 
shown in Fig. [TJ with th e most prominent H II regions labeled, following the numbering 
system of iPustilnik et al. I (120051 ). It is seen from the g — i image that the light of DDO 68 is 
dominated by blue young stellar populations (in black) without evidence for a red halo. This 
suggests that DDO 68 m ay be a young system, with an age less than 1 Gyr, as discussed by 



Pustilnik et al. I (j2008al ). The several compact red objects in the image (in white) are either 



foreground Galactic stars or distant background galaxies. 



2.2. Spectroscopy 

New high signal-to-noise ratio optical spectra of DDO 68 were obtained using the 3.5 m 
Apache Point Observatory (APO)§ telescope and the 6.5 m MMT0 on the nights of 2 February 
2008 and 28 March 2008, respectively. The APO observations of DDO 68 were made with the 
Dual Imaging Spectrograph (DIS), which covers both the blue and red wavelength ranges. A 
l'.'5x360" slit was used. In the blue range, we use the B400 grating with a linear dispersion 
of 1.83 A/pix and a central wavelength of 4400A, while in the red range we use the R300 
grating with a linear dispersion of 2.31 A/pix and a central wavelength of 7500A. The above 
instrumental set-up gave a spatial scale along the slit of 0'.'4 pixel -1 , a spectral range of 
~3600 - 9600A and a spectral resolution of 7A (FWHM). The slit was oriented in such a 
way as to include both regions 3 and 4 (hereafter DDO 68-3 and DDO 68-4) (Figured]). The 
MMT observations of DDO 68 were made with the Blue Channel spectrograph. We used a 
l'.'5xl80" slit and a 800 grooves/mm grating in first order. The above instrumental set-up 
gave a spatial scale along the slit of 0'.'6 pixel -1 , a scale perpendicular to the slit of 0.75A 
pixel -1 , a spectral range 3200 - 5000A and a spectral resolution of 3A (FWHM). We show 



2 IRAF is distributed by National Optical Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 

3 The Apache Point Observatory 3.5-meter telescope is owned and operated by the Astrophysical Research 
Consortium. 

4 The MMT is operated by the MMT Observatory (MMTO), a joint venture of the Smithsonian Institution 
and the University of Arizona. 
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the MMT slit location superposed on an expanded image of DDO 68-3 in Fig. [2h. The 
star-forming region is nearly unresolved. However, the brightness distribution along the slit 
suggests that the LBV star is contained within it (Fig. [2b). 

As for PHL 293B, we did not obtain new observations but used a high spectral resolution 
spectrum of PHL 293B in the archives of the European Southern Observatory (ESO). The 
spectrum was obtained with the UVES spectrograph on the 8m VLT@ during the night of 8 
November 2002, in both the blue arm (grating CD^2, central wavelength 3900A, slit l"x8") 
and the red arm modes (grating CD^3, central wavelength 5800A, slit l"xl2"), giving a 
wavelength range 3100 - 6800A. The resolving power was ~ 80,000, resulting in a spectral 
resolution ~ 0.2A. The VLT slit superposed on a g band SDSS image of PHL 293B is shown 
in Fig. [2b. As in the case of DDO 68-3, the brightness distribution along the slit (Fig. 
[2bl) suggests that the LBV star resides within the bright, compact and barely resolved star- 
forming region. The ESO spectrum of PHL 293B was supplemented by the SDSS spectrum, 
obtained on 22 August 2001, with a wavelength range 3800 - 9200 A and a spectral resolution 
of ~ 3A. The details of the spectroscopic observations are also given in Table [TJ 

For both APO and MMT observations, several Kitt Peak and CTIO IRS spectroscopic 
standard stars were observed for flux calibration. Spectra of He-Ne-Ar comparison arcs were 
obtained at the beginning of each night for wavelength calibration. As for the high spectral 
resolution VLT/UVES observations, calibration was done by observing the standard star 
Feige 110, and wavelength calibration was done using the spectrum of a Thorium comparison 
arc. 

The data reduction pr o cedur es of the spectroscopic observations are the same as de- 



scribed in llzotov fc Thuan I (120071 ). Briefly, the two-dimensional spectra were first bias sub- 
tracted and flat-field corrected with IRAF. We then use the IRAF software routines IDEN- 
TIFY, REIDENTIFY, FITCOORD, TRANSFORM to perform wavelength calibration and 
correct for distortion and tilt for each frame. Background subtraction was performed us- 
ing the routine BACKGROUND. The level of background emission was determined from the 
closest regions to the galaxy that are free of galaxian stellar and nebular line emission, as well 
as of emission from foreground and background sources. A one-dimensional spectrum was 
then extracted from the two-dimensional frame using the APALL routine. Before extraction, 
separate two-dimensional spectra of each object were carefully aligned using the location of 
the brightest part in each spectrum, so that spectra were extracted at the same positions 
in all subexposures. We then summed the individual spectra from each subexposure after 



5 Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme 
ID 70.B-0717(A) 
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manual removal of the cosmic rays hits. 

The resulting MMT spectra of DDO 68-3 and DDO 68-4 are shown in Fig. [3J The 
spectrum of DDO 68-3 is strikingly different from normal spectra of H II regions: it shows 
strong broad hydrogen and helium emission lines with P Cygni profiles superposed on narrow 
nebular emission lines (Fig. [3^). No emission lines of heavy elements (e.g., Fe II and N ill 
lines) were detected, in contrast to high-metallicity quiescent LBV spectra which show the 
presence of many such lines. The absence of permitted heavy element emission lines in the 
visible spectrum of the LBV star in DDO 68-3 is most likely due to its very low metallicity, 
but probably also to temperature effects in the pseudo-photosphere of the LBV star as it 
is likely in an eruptive phase. No broad emission is detected in DDO 68-4 (Fig. [lb). The 
[O ill] A4363 emission line is detected in both DDO 68-3 and DDO 68-4, allowing a direct 
determination of the electron temperature and element abundances. 

The high-resolution VLT spectrum of PHL 293B is shown in Fig. [H Broad emission 
and P Cygni profiles are clearly seen for the Ha, H/3 and H7 lines, superposed on narrow 
nebular emission lines. The broad emission is not evident in the higher-order hydrogen lines, 
even if blue absorption is present (marked by arrows in Fig. H]) . However, at variance with 
DDO 68, no broad He I emission is seen in the spectrum of PHL 293B. Similar to DDO 
68, no permitted heavy element emission lines with P Cygni profiles were detected in the 
spectrum of PHL 293B, again a consequence of the low metallicity of the LBV star. 



3. RESULTS AND DISCUSSION 
3.1. Photometry of the LBV star in DDO 68 

We now examine the temporal variation of the brightness of the LBV star in DDO 68, 



i n an a ttempt to narrow down the time of its outburst. It was first seen by lPustilnik et al. 



( I2008bl ) thanks to the presence of broad emission in the spectrum of DDO 6 8-3 obtained on 



11 Jan uary 2008. Earlier spectroscopic observations of the same region by iPustilnik et al 



( 120051 ) in 2005 did not reveal anything abnormal. This implies that outburst occurred 
sometime during the 2005 - 2007 period. We can use the SDSS and our 2.1m KPNO imaging 
observations of DDO 68 to narrow down this range further. With the IRAF APPHOT 
package, we have measured on both sets of images g and i magnitudes of the H II regions in 
DDO 68 that are labeled in Fig. Qja. The SDSS images were obtained on 16 April 2004 and 
the KPNO images on 9 February 2007. The results of our photometric measurements are 
shown in Table [2J Examination of this table shows that the g magnitudes of regions 1-6 
derived from KPNO and SDSS images are in very good agreement when measurement errors 
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are taken into account. Thus, the g magnitude difference between the two sets of images for 
DDO 68-3 (where the LBV star has been identified spectroscopically) is only -0.06 mag and 
can be attributed to photometric errors. The agreement between the i magnitudes is not as 
good, primarily because of the noisy short exposure SDSS % image. We thus conclude that 
no drastic luminosity change has occurred in DDO 68-3 up to 9 February 2007, and that the 
LBV outburst has occurred sometime between 9 February 2007 and 11 January 2008. 

As for PHL 293B, except for the SDSS images, we do not possess other images of it 
taken at subsequent times. All we can say is that the LBV outburst in it occurred prior to 
22 August 2001, the date of the SDSS spectrum, as it already showed broad emission by 
then. Clearly, imaging of both DDO 68 and PHL 293B is urgently needed to constrain the 
magnitude jump of the LBV outburst in the two BCDs. Are both eruptions "normal" LBV 
outbursts with 5m ~ 1-2 mag, or major outbursts with 5m > 3 mag? 



3.2. Element abundances 



We now derive element abundances for both objects from the narrow emission line fluxes. 
These fluxes have been measured using Gaussian fitting with the IRAF SPLOT routine. They 
have been corrected for both extinction, using the reddening curve of IWhitford I (119581 ). and 
underlying h ydrogen stella r abso rption, derived simultaneously by an iterative procedure as 
described by llzotov et al. I (119941 ) and using the observed decrements of the narrow hydrogen 
Balmer lines. It is assumed in this procedure that hydrogen line emission is produced only 
by spontaneous transitions in recombination cascades, i.e. we neglect possible collisional 
excitation. Such a situation usually holds in low-density H II regions ionized by stellar 
radiation such as those considered here. The extinction-corrected fluxes 100xI(A)//(H/?) 
of the narrow lines for each galaxy, together with the extinction coefficient C(H/3), the 
equivalent width of the FL3 emission line EW(H/3), the H{3 observed flux F(H./3) and the 
equivalent widths of the underlying hydrogen absorption lines EW(abs) are given in Table [3j 
The physical conditions and element abundances of DDO 68-3 (APO and MMT observations) 
and DDO 68-4 (MMT observations) and of the H II region in PHL 293 B (VLT observations) 
are derived from the narrow line fluxes following llzotov et al. I (120061 ) . Briefly, we adopt a 
two- zone photoionized H II region model: a high-ionization zone with electron temperature 
T e ([0 m]), where O ill, Ne ill, and Ar iv originate, and a low- ionization zone with electron 
temperature T e ([0 n]), where O II, N II, S II and Fe ill originate. The temperature T e ([0 
ill]) is derived from the [O ill] A4363/(A4959+A5007) flux ratio. This temperature is used 
for the determination of the 2+ and Ne 2+ ionic abundances. The tem perature T e (\Q n]) is 
obtained from the relation between T e ([0 n] and T e ([0 ill] derived by llzotov et al. I (120061 ) 
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by fitting photoionization models. This temperature is used for the determination of the 
+ and N + ionic abundances. The electron number density N e is derived from the [S n] 
A6717/6731 flux ratio when it is available, and is set to 10 cm -3 when it is not. The total 



eleme nt abundances are obtained using ionization correction factors given by llzotov et al. 
(120061 ) . The element abundances are given in Table HI 



Within the errors, the oxygen abundance derived for DDO 68-3 from the APO spectrum 
is in agreement with the one derived from the MMT spectrum. Note, however, that the 
abundance determination from the lower signal-to-noise ratio and lower resolution APO 
spectrum is much more uncertain. The derived values 12 + log O/H = 7.15 ± 0.04 (DDO 
68-3) and 12 + log O/H = 7.16 ± 0.09 (DDO 68-4) are in very goo d agreement with 
12 + log O/H = 7.14 ± 0.03, 7.13 ± 0.07 and 7.21 ± 0.14 derived bv llzotov fc Thuan 
(120071 ) for regions 1, 2 and 6, respectively. T hey are also con s istent wit h the abundance 



determ inations for other regions in DDO 68 by lPustilnik et al. I (120051 ) and iPustilnik et al. 



(j2008bl ) from lower signal-to-noise ratio spectra. There appears not to be a large abundance 
variation among various H II regions in DDO 68 over a spatial scale of ~ 1 kpc. The relative 
abundance homogeneity puts stringent constraints on the physics of metal diffusion and 
mixing pr ocesses in this object. T he overall conclusion is that we confirm and reenforce the 
finding by llzotov fc Thuan I (120071 ) that DDO 68 is the second lowest-metallicity emission-line 
galax y known, after SBS 0335-052W which has 12 + log O/H = 7.12 ± 0.03 (llzotov et al 



20051 ). This makes the discovery of a LBV star in DDO 68 particularly important. Being the 
lowest metallicity LBV known, its study will allow to put stringent constraints on evolution 
models of metal-deficient massive stars. 

The oxygen abundance of PHL 293B derived from the VLT spectrum is significantly 
higher than those of the H II regions in DDO 68, being equal to 12 + log O/H = 7.72 ± 0.01. 
But if DDO 68 is excluded, it is still considerably lower than the oxygen abundances of all 
other galaxies with det ected LBVs. Within the errors, it is consistent with 12 + log O/H = 
7.66 ± 0.04 derived by llzotov et al. I (120071 ). from the lower resolution SDSS spectrum. It 
is also interesting to note that the nitrogen-to-oxygen abundance ratio log N/O = -1.55 ± 



0.02 is typical of low-metallicity emission-line galaxies (llzotov et al. 1120061 ) 



3.3. Broad emission 

Using Gaussian fitting, we have measured the fluxes of the broad components of the 
Balmer H lines in the spectra of DDO 68-3 and PHL 293B, and of the He I lines in the 
spectrum of DDO 68-3 after subtraction of the narrow component. The extinction-corrected 
fluxes of the broad Ha, H/3 and H7 lines together with the fluxes of the narrow [O ill] A4959 
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and 5007 nebular lines, derived from the MMT observations of DDO 68-3 and from the VLT 
observations of PHL 293B, are shown in Table O We also present in the same table the 
fluxes of the same lines, derived from the APO spectrum for DDO 68-3 and from the SDSS 
spectrum for PHL 293B. Comparison of the data for DDO 68-3 obtained at different epochs 
(2 February 2008 for the APO spectrum and 28 March 2008 for the MMT spectrum) reveals a 
flux increase of a factor of ~ 2 for H/3 over the course of nearly 2 months. This is likely a real 
change as the fluxes of the nebular emission lines [O ill] A4959, 5007 remain approximately 
constan t over this period. This conclusion is further supported by the observations of DDO 



68-3 bv lPustilnik et al. I (I2008H ). They derived a broad H/3 flux of 5.46 xlO -16 erg s _1 cm" 2 



on 11 January 2008, or ~ 2 times lower than the broad H/3 flux observed on 2 February 2008 
at the APO and ~ 4 times lower than the one observed on 28 March 2008 at the MMT, 2.5 
months later. It thus appears that the LBV star in DDO 68 is experiencing a significant 
increase of emission line brightness. In Fig. [SJ we compare the MMT spectrum (thick solid 
line) with t he APO spectrum (thin solid line). The continuum level before the outburst as 
observed bv lPustilnik et al. I (j2008bl ) (their Fig. 1, middle panel), is shown by a dotted line. 
It can be seen that the line flux brightening is accompanied by an increase in the continuum 
level. The continuum level of the MMT spectrum at A4686, the av erage wavelength of the q 



band, is higher by a factor of ~ 3 than the continuum level of the iPustilnik et al. I (2008bJ ) 



spectrum, taken before the outburst. Adopting an apparent magnitude g = 20.17 for DDO 
68-3 before the erupti on (Table [2]) and a dista nce modulus of 30 mag (which corresponds to 
a distance of 10 Mpc; IPustilnik et al. Il2008bl ). the g absolute brightness of the LBV star in 
DDO 68-3 is ~ —11 mag on 28 March 2008. If we adopt the smaller distance of 6.24 Mpc given 
in the NEL@, then the g absolute magnitude of the LBV star is ~ —10 mag. In either case , 



the a bsolute bri ghtness is in the r ange o f those of known LBVs ([Humphreys fc Davidson 



19941 ) . Note that IPustilnik et al. I (j2008bl ) also found an increase in the continuum level of a 
factor of ~ 2 between 11 January 2008 and 4 February 2008 (their Fig. 2). We plan to obtain 
both images and spectra of DDO 68 at regular intervals to monitor the LBVs brightness 
evolution in time. 

As for PHL 293B, the fluxes of the broad H lines appear also at first glance to vary 
with time: they decrease by a factor of nearly 2 between the SDSS (22 August 2001) and 
VLT (8 November 2008) observations. However, a more detailed look casts doubt on that 
interpretation: the fluxes of the nebular [O ill] A4959, 5007 emission lines also decrease by 
a similar factor of ~ 1.5. Therefore, these variations may not be real and could be caused 
by aperture differences between the VLT and SDSS observations. A narrow slit of 1" width 



6 NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under contract with the National Aeronautics and Space Administration. 
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was used in the VLT observations while a round 3" aperture was used during the SDSS 
observations. In contrast to the VLT slit, the larger SDSS aperture probably included all 
the light from the LBV star. 

The Ha luminosities of the LBV stars (Table [5]) have been calculated adopting distances 
from the NED D = 6.24 Mpc and 22.7 Mpc for DDO 68 and PHL 293B, respectively. These 
distances have been obtained from the radial velocities corrected for Virgo Infall with a 
Hubble constant of 73 km s -1 Mpc -1 . The derived luminosity L(Ha) = 9.4xl0 36 erg s _1 
(Table [5]) for the LBV in DDO 68 from the APO spectrum is significantly lower than L(Ha) 
= 3.4 x 10 38 erg s _1 of the LBV in PHL 293B as derived from the SDSS spectrum. We have no 
data for the Ha flux for the DDO 68 star from the MMT observations as they did not cover 
the red wavelength range. However, if we accept that the Ha flux has brightened by a factor 
of ~ 2 in the period between the APO and MMT oservations as the H/3 flux and that the 



distance to DDO 68 can be as high as 10 Mpc (see discussion in iPustilnik et al. Il2008bl ). we 
conclude that the Ha luminosity of the LBV in DDO 68 can be as high as 5xl0 37 erg s -1 on 
28 March 2008. Then the Ha luminosities of the LBVs in DDO 68 and PHL 293B compare 
well with L( H.a) ~ 10 38 erg s -1 of the VI star in NGC 2363, the third lowest-metallicity 



LBV known flDrissen et al. Il2001[ ) after those in PHL 293B and DDO 68. 



3.4. Stellar winds 



In Table [6], we show the terminal velocities i>oo of the stellar winds associated with the 
outbursts and the FWHMs of the broad H and He I lines of the DDO 68 and PHL 293B 
LBVs, as derived from the high-resolution MMT and VLT spectra. The stellar wind terminal 
velocity is derived from the wavelength difference between the blue absorption minimum 
and the broad emission line maximum. For both LBVs, the i>oo's and FWHMs derived from 
different lines are in very good agreement. Our derived FWHMs for the LBV in DDO 68- 3 
of ~ 1000 km s -1 are in good agreement with those derived by IPustilnik et al. I (l2008bl ). 



However, while our FW HMs are consistent from line to line (Table E]), those derived by 
Pustilnik et al. I (l2008bl ) show a large spread, mainly because of their lower signal-to-noise 
ratio spectra. The FWHMs measured for the LBV in PHL 293B are ~ 600 km s _1 , lower 
than those of the LBV in DDO 68 by a factor of ~ 1.7. However, the FWHMs of both 
LBVs are significantly broader than t hose observed in the s pectra of high-metallicity LBVs 



which are about 100-200 km s 1 (e.g. iLeitherer et al. II19941 ). These differences are probably 



not all due to metallicity effects. Part of the FWHM differences may be explained by the 
fact that the high-metallicity LBVs are observed during their quiescent phase, while our 
low- metallicity LBVs are likely observed during their eruptive phase. 
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The same high velocities are found when Voo is considered. Indeed, one of the most 
striking features of low-metallicity LBVs is that the terminal velocities of their winds are 
significantly higher than those of higher metallicity LBVs. The range of terminal velocities 
in higher metallicity LBVs, including VI in NGC 2363, is 2 00 - 400 km s -1 . For example, 



the galactic LBV P Cygni has a of only 185 km s -1 (ICrowther 1 120041 ). while in the 
low-metallicity LBVs is ~ 800 km s _1 . Probably, selection effects may play a role here: 
generally, emission-line galaxies are observed spectroscopically at low resolution, so that 
their spectra do not allow the finding of LBVs with low v^. In particular, with the spectral 
resolution of ~ 2-3 A of the SDSS spectra, it is unlikely that we can use them to find LBV 
with terminal velocities < 200 km s _1 . On the other hand, no high- metallicity LBV was 
ever found to have high terminal velocities. Most likely, metallicity does play a role and 
the stellar winds in low-metallicity LBVs are driven by a mechanism which is different from 
the one operating in high-metallicit y LBVs, thought t o be th e radiation pressure due to the 



permitted lines of heavy elements. ISmith fc Owocki I (120061 ) have considered a continuum 



driven mechanism which could be more efficient at low metallicities than the line-driven 
mechanism. It is likely that winds in low-metallicity LBVs are more transparent and less 
dense, which results in higher terminal velocities. 

New high spatial resolution spectral and photometric observations are necessary to shed 
light on the physical properties and evolution of low-metallicity LBVs. They will strongly 
constrain evolutionary models of low-metallicity massive stars on the post-main-sequence 
stages of their evolution. This will in turn increase our understanding of the evolution of the 
first generation of massive stars, when the gas was pristine. 



4. CONCLUSIONS 

We have studied here the broad line emission of two luminous blue variable (LBV) stars 
discovered in two low-metallicity blue compact dwarf (BCD) galaxies, DDO 68 and PHL 
293B. We have arrived at the following conclusions: 

1. The oxygen abundances in region 3 of DDO 68 (DDO 68-3) and in PHL 293B are 
respectively 12+log O/H = 7.15 ± 0.04 and 12+log O/H = 7.72 ± 0.01. These two BCDs 
are thus the lowest-metallicity galaxies with detected LBV stars. PHL 293B is also the most 
distant galaxy where a LBV star is seen. 

2. Photometric observations of DDO 68-3 show that the outburst in its LBV star has 
occurred sometime between 9 February 2007 and 11 January 2008. 

3. Broad H and He I emission lines with P Cygni profiles are seen in the spectrum of 
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the LBV star in DDO 68. On the other hand, only H broad emission lines are detected 
in the spectrum of the LBV in PHL 293B. In both LBVs, no heavy element emission line 
such as Fe II was detected, presumably because of their low metallicities. The broad Ha 
lumin osities in both LBVs compare well with the on e in the higher-metallicity LBV in NGC 
2363 JPrissen et al. I ll997L I2OO1I : IPetit et al. Il2006h . We find a strong increase of the Ha 
luminosity of the LBV in DDO 68: during the period from 11 January 2008 to 28 March 
2008, it brightened by a factor of ~ 4. 

4. One of the most striking features of both low-metallicity LBVs is the high terminal 
velocities Voo of their stellar winds. Their Voo are ~ 800 km s _1 , several times greater than 
the of ~ 100 - 200 km s -1 in their high-metallicity counterparts. Probably, metallicity 
plays a role and winds in low-metallicity LBVs are driven by a mechanism which differs from 
the one operating in high-metallicity LBVs. 
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Table 1. Journal of Observations 



Telescope 


Object 


Date 


Wavelength range 


Slit 


Exposure 


Airmass 


2.1m KPNO 


DDO 68 


09 Feb 2007 


g-band 




1800 


1.01 


2.1m KPNO 


DDO 68 


09 Feb 2007 


i-band 




1800 


1.05 


2.5m Sloan 


DDO 68 


16 Apr 2004 


g-band 




54 




2.5m Sloan 


DDO 68 


16 Apr 2004 


i-band 




54 




3.5m APO 


DDO 68 


02 Feb 2008 


3600-9000 


l'/5x300" 


2700 


1.70 


6.5m MMT 


DDO 68 


28 Mar 2008 


3200-5200 


l'.'5xl80" 


4500 


1.08 


2.5m Sloan 


PHL 293B 


22 Aug 2001 


3800-9200 


3" (round) 


2400 


1.25 


8m VLT 


PHL 293B 


08 Nov 2002 


3100-4500 


l"x8"(blue) 


3300 


1.10 








4800-6800 


l"xl2"(red) 


3300 


1.10 
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Table 2. 


Photometry of DDO 68 






Region 


2.1m KPNO 


2.5m Sloan 




Ai a 


9 i 


9 i 




09 Feb 2007 


16 Apr 2004 






1 


19.44±0.02 20.27±0.05 


19.41±0.04 20.46±0.18 


+0.03±0.04 


-0.19±0.19 


2 


20.44±0.03 21.38±0.09 


20.40±0.07 21.05±0.31 


+0.04±0.08 


+0.33±0.32 


3 


20.17±0.02 20.63±0.06 


20.23±0.05 21.03±0.29 


-0.06±0.05 


-0.40±0.30 


4 


20.19±0.04 21.04±0.10 


20.18±0.07 21.66±0.54 


+0.01±0.08 


-0.62±0.55 


5 


19.13±0.01 19.41±0.02 


19.08±0.02 19.54±0.08 


+0.05±0.02 


-0.13±0.08 


6 


20.40±0.04 20.84±0.06 


20.38±0.07 20.40±0.15 


+0.02±0.08 


+0.44±0.16 



a magnitude difference between KPNO and SDSS measurements 



Table 3. Intensities of narrow emission lines 



Ion 100x/(A)//(H/?) 100x/(A)//(H/3) 100 x /(A) //(H/?) lOOx J(A)/J(H/3) 



Galaxy 





DDO 68-3 


DDO 68-3 


DDO 68- 


4 


PHL 293B 




(MMT) 


(APO) 


(MMT) 




(VLT) 


3727 [0 n] 

L J 


50.62 ± 1.04 


57.83 


± 2.59 


94.10 ± 


1.70 


50.21 ± 0.84 


3868 [Ne ill] 


10.07 ± 0.40 


6.62 


± 0.89 


5.22 ± 


0.09 


48.49 ± 0.77 


3889 He I + H8 


17.74 ± 0.62 


19.81 


± 1.91 


18.13 ± 


1.20 


21.27 ± 0.41 


3968 [Ne ml + H7 


16.35 ± 0.59 


21.39 


± 1.78 


15.72 ± 


1.13 


31.70 ± 0.54 


4101 E5 


26.73 ± 0.68 


28.39 


± 1.67 


25.56 ± 


1.14 


27.00 ± 0.46 


4340 H 7 


49.34 ± 0.98 


46.39 


± 1.70 


47.71 ± 


1.27 


48.25 ± 0.74 


4363 [O ill] 


2.82 ± 0.25 


2.28 


± 0.53 


1.64 ± 


0.38 


14.55 ± 0.27 


4471 He I 


3.28 ± 0.25 










3.75 ± 0.10 


4861 HP 


100.00 ± 1.70 


100.00 


± 2.30 


100.00 ± 


1.95 


100.00 ± 1.44 


4959 [0 ill] 


39.56 ± 0.71 


27.37 


± 0.83 


23.79 ± 


0.72 


195.13 ± 2.79 


5007 [0 ill] 


115.29 ± 1.84 


92.85 


± 2.03 


68.36 ± 


1.40 


584.35 ± 8.33 


5876 He I 












10.07 ± 0.17 


6563 Ha 




276.69 


± 5.77 






277.28 ± 4.31 


6583 [N II] 












1.82 ± 0.06 


6678 He I 












2.76 ± 0.05 


6717 [S n] 




4.31 


± 0.30 






5.04 ± 0.11 


6731 [S n] 




3.00 


± 0.27 






3.68 ± 0.08 


C(H0) 


0.000 




0.005 


0.000 


0.080 


EW(H/?) A 


16 




85 




36 


37 


F(H/3) a 


21.9 




29.5 




9.2 


71.5 


EW(abs) A 


0.00 




1.55 




0.80 


0.05 



a in units 10 16 erg s 1 cm 2 . 
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Table 4. Element abundances 



Object 12+logO/H log N/0 log Ne/O 



DDO 68-3 (MMT) 


7.15 


± 


0.04 


-0.81 


± 


0.07 


DDO 68-3 (APO) 


7.08 


± 


0.09 


-0.93 


± 


0.22 


DDO 68-4 (MMT) 


7.16 


± 


0.09 


-1.02 


± 


0.25 


PHL 293B (VLT) 


7.72 


± 


0.01 


-1.55 ± 0.02 -0.74 


± 


0.02 
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Table 5. Fluxes and luminosities of the broad emission lines at different epochs 



Telescope 


Date 


/(H 7 ) a 


J(H/3) a J(H«) a 


L(Ha) b 


J(4959) a 


J(5007) a 


DDO 68-3 


APO 
MMT 


2008-02-07 
2008-03-28 


10.9±0.4 


9.7±0.5 20.2±0.4 
19.7±0.3 


9.4xl0 36 


8.4±0.3 
8.8±0.2 


26.6±0.4 
25.5±0.4 


PHL 293B 


SDSS 
VLT 


2001- 08-22 

2002- 11-08 


7.5±0.2 


26.6±4.5 97.8±6.3 
16.7±0.3 55.2±0.6 


6.0xl0 38 
3.4xl0 38 


196.5±6.8 
139.9±1.5 


603.8±20.1 
412.8±4.0 



a In units 10 16 erg s 1 cm 2 . Flux errors are derived taking into account photon statistics in 
non-flux calibrated spectra. 

b In erg s _1 . 
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Table 6. Terminal velocities and FWHMs of the broad emission lines in DDO 68-3 and 

PHL 293B 





DDO 6! 


3-3 (MMT) 


PHL 293B (VLT) 


Line 


Voo a 


FWHM k a 




FWHM fer a 


6563 Ha 






859 


719±10 


5016 He I 


826 


974±62 






4861 H/3 


801 


1008±10 


819 


649±12 


4713 He I 


807 


1004±89 






4471 He I 


803 


1053±54 






4340 H 7 


835 


1122±20 


871 


681±14 


4101 m 


788 


1137±33 


735 


474±30 


3968 H7 


762 


1100±40 






3888 H8 


732 


1015±31 







a In km s 1 . 
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Fig. 1. — a) 2f5 x2.5 2. 1m KPNO (7 image of DDO 68. The H II regions are labeled following 
Pnstilnik et XI (hoOfih . b) 2.1m KPNO 5 - i image of DDO 68. 
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-5 5 -5 5 

Distance along the slit (arcsec) Distance along the slit (arcsec) 



Fig. 2.— a) MMT slit superposed on a KPNO 2.1m g band image of DDO 68-3. b) VLT slit 
superposed on a SDSS g band image of PHL 293B. c) Normalized B.(3 brightness distribution 
(solid line) and normalized continuum distribution (dashed line) of DDO 68-3 along the MMT 
slit, d) The same brightness distributions as in c) but in PHL 293B along the VLT slit. 
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6 - 



o L 



DDO 68-3 





a) 



3500 



4000 



4500 
Wavelength (A) 



5000 



DDO 68-4 



2 - 



1 - 




b) 



3500 



4000 



4500 
Wavelength (A) 



5000 



Fig. 3.— Redshift-corrected 6.5 MMT spectra of DDO 68-3 (a) DDO 68-4 (b). The emission 
lines are labeled in panel a). 
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Fig. 4. — Redshift-corrected archival 8m VLT/UVES spectrum of PHL 293B. The locations 
of blue-shifted absorption in lines with P Cygni profiles are shown by arrows. 
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4500 
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Fig. 5. — Comparison of the MMT spectrum (thick solid line) and of the APO apectrum 
(thin solid line) of DDO 68-3 shows an increase of the continuum level of a factor of ~ 3. 
Th e dotted line shows the continuum level before the LBV outburst (middle panel of Fig.l 
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